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VACANCY FORMATION IN GOLD UNDER HIGH PRESSURE 
by Hubert H. Grimes 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

High-purity gold wires were resistively heated and quenched in a 
pure argon gas environment at pressures up to 6 kilobars. Quench rates 
of 10^ to 5x10^ °C/sec were obtained by abruptly reducing the heating 
current. The resistance of both the specimen and an unquenched dummy 
specimen were measured in situ at -38° C with conventional circuits. 

The logarithm of the quenched-in resistance decreased linearly with in- 
creasing pressure at a constant temperature of 600° C. From these data, 
a formation volume of 0.52±0.07 atomic volume was derived, which was as- 
sumed to be associated with the formation of single vacancies. For at- 
mospheric pressure quenches from below 725° C, Arrhenius plots of the 
quenched-in resistivity for various quenching temperatures yielded an 
average formation energy for single vacancies of 0.93 ev. Above 725° C, 
the quenched-in resistivity drops below this line consistent with the 
complex annealing behavior during the quench observed by others for 
higher temperature quenches. The measured formation volume indicated 
a degree of relaxation around a vacancy that is shown to be in accord 
with previously calculated lattice relaxations obtained 
function theory. 

INTRODUCTION 

Large elastic strains have been shown to have a significant effect 

O £ 

on self -diffusion in several metals. 




Compressive strains resulting 


2 


* from the application of hydrostatic pressures of 5 to 10 kilobars have 
reduced the diffusion coefficient by as much as an order of magnitude. 
Theory predicts that the principal manifestation of the strain arises 
in the exponential part of the diffusion equation associated with the 
energy required to form a vacancy in the deformed crystal and with the 
energy an atom requires to surmount a potential barrier along the diffu- 
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sion path. But this theory does not enable one to predict a priori the 
relative importance of these two effects. However, several quenching 
and annealing experiments on the noble metals have resulted in the inde- 
pendent determination of the formation energy and migration energy of 
vacancies in unstrained materials. These techniques have been extended 

O 

to strained systems for vacancy annealing in gold, and at the time of 

Q 

this writing, for quenching in gold. Other methods, the change in mac- 
roscopic length^ and calorimetric measurement of the energy stored in 
a quenched crystal, have been used to determine, indirectly, the ef- 
fect of strain on vacancy formation. 

The technique reported herein is the most direct; it involves de- 
termination of the vacancy concentration by a measurement of the excess 
electrical resistivity quenched into gold wires that are subjected to 
large hydrostatic pressures. The data permit the calculation of not 
only the energy of formation of the vacancy but also the volume change 
of the crystal accompanying the formation of a vacancy. The volume ob- 
tained indicates the degree of relaxation of the lattice atoms around 
the vacancy. This relaxation is compared with calculated estimates of 
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this quantity obtained using simple pairwise potential function theory. 
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THEORY 


It has been shown that by rapidly quenching fine metal wires from 
elevated temperature to room temperature or below, a large fraction of 
the point defects in equilibrium with the crystal lattice at the high 
temperature, can be retained. 10 The measure of this defect supersatura- 
tion at the low temperature is reflected in an increase in the electri- 
cal resistivity of the wire. Annealing of these defects in gold has 
indicated that defects quenched from below 700° C, at quench rates of 
10 3 to 10 5 °C/sec or more, are predominately single vacancies. Above 
this temperature, more complex defect annealing occurs and has been the 
subject of several theoretical and experimental studies. 13-13 For sim- 
plicity of interpretation, therefore, the quenches under high pressures 
in this study were done at 600° C. 

The number of vacancies n^, which are in equilibrium with a crys- 
tal lattice at absolute temperature T, can be expressed as 

-AG f /kT 


r\y = Ne 


( 1 ) 


where N is the total number of possible lattice sites at which a va- 
cancy can exist and AGf is the change in the Gibbs free energy of the 
crystal accompanying the formation of a vacancy. Using the thermody- 
namic relation that the Gibbs free energy is related to the formation 
energy AH f , or enthalpy, and to the entropy of formation AS f by the 
expression 

£G f = AH f - T AS f (2) 


and also using the observation that the fractional vacancy concentration 
quenched into a crystal from a temperature Tq at pressure p is pro- 
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• portional to the resistivity increase £p v on quenching, one obtains 
for the quenched-in resistivity, 



-Z£„/kT n AS„/k -AH„/kT n 
Ae r * * Ae f e f Q 


(3) 


where A is a constant independent of temperature and pressure. It 
follows from Eq. (3) that In Ap y is linear in l/T q with a slope 

6 (In Ap y ) Ni f £H f 
6(1/T q ) = " “ = ' “k~ 

if AS|. is also assumed to he independent of temperature. This last 
assumption has been justified experimentally, and the energy of forma- 
tion of a vacancy may be derived from the slope of such a curve. 

At constant temperature the derivative of the Gibbs free energy 
with pressure is a volume, which for dGf has the interpretation of 
the volume change a crystal experiences on formation of a vacancy. 

This derivative, from use of Eq. (3), is 



Accordingly, the formation volume can be determined from a plot of 
In Ap v against pressure for quenches from a constant temperature. 

EXPERIMENTAL 

The large hydrostatic pressures required in this experiment were 
obtained with a Harwood Engineering Company three-stage gas pumping sys- 
tem capable of attaining pressures up to 13 kilobars in a cylindrical 
vessel with internal dimensions of 2 inches in diameter and 8 inches 
high. , High-purity argon gas (99.995 percent) was used as the high- 
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pressure fluid because it is chemically inert and has favorable thermo- 
dynamic properties. Gas chromatographic analysis of the argon confirmed 
the listed purity. The pressure system was operated with a minimum of 
lubrication to avoid thermal decomposition of oil on the specimen at 
elevated temperatures. The high pressures were measured with a cali- 
brated Manganin gage , linear to within 1 percent in the pressure range 
used. 

The specimen holder used for both the high- and low-pressure 
quenches is shown in Fig. 1. It is essentially an aluminum housing that 
encloses horizontally suspended specimen and dummy specimen wires. This 
enclosure is supported in good thermal contact with a 7/8 in. (O.D.), 
l/4 in. (I.D. ) steel tube that is closed on the lower end. Ethanol 
cooled to dry-ice - ethanol bath temperature was circulated through this 
steel "cold thimble" through a specially designed pressure vessel clo- 
sure. In this way it was possible to cool the specimen enclosure even 
when the vessel was at high pressure. 

The specimen enclosure was built up in sections to facilitate the 
insertion and suspension of the specimen and dummy specimen and the 
mounting of thermocouples. In the finished construction the specimen 
and dummy specimen wires remained suspended at the center of individual 
1/4 -in. cross-sectional annular cavities surrounding the cold thimble. 
With this arrangement, the temperature of the wires could be precisely 
maintained near -38° C by heat exchange with the cold walls of the 
cavities. The enclosure was provided with small baffled vent holes to 
permit passage of the high pressure gas. The entire specimen holder 
assembly was surrounded with a Lucite jacket that both insulated the 
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* cold assembly from the pressure vessel wall and reduced the dead space 
in the vessel. 

A conventional- circuit, D shown schematically in Fig. 2, was used 
in the measurement of the quenched-in resistance. A differential meas- 
uring technique, used to improve sensitivity, compared the resistance of 
the quenched specimen to an unquenched dummy specimen. The dummy speci- 
men also served to compensate for ambient temperature changes of the as- 
sembly before and after the quench. 

To make a measurement, the bridge circuit, including the specimen, 
Rg the dummy, R^ and two 1-ohm standard resistors R was balanced 
by adjusting the variable resistors R A and Rg. In the balanced 
bridge, the measuring currents in the specimen and the dummy could be 
maintained equal to 1 part in 10 5 . The resistance difference between 
the specimen and dummy specimen was then determined by measurement of 
the potential differences E]_ and Eg with a modified White potentiom- 
eter accurate to 10"® volt . Errors due to constant thermal emf ' s were 
eliminated by reversing the measuring current flow in the circuit. The 
resistance difference was then given by 


E, - E, 


R g - R D = 


where E-j_ and Eg are averages of the forward and reverse readings. 
The measuring current I (nominally 10 ma in each arm of the bridge) 
was determined by measuring the potential drop across one of the 1-ohm 
standard resistors. This value of measuring current was found to mini- 
mize heating of the specimen and the dummy specimen, while permitting 
the quenched-in resistance to be determined to within 1 percent of the 
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* total value for most quenches. 

The specimens and dummy specimens were, in all cases, fabricated 
from 0.003-in. -diameter gold wire from the same lot, quoted by the sup- 
plier 17 as better than 99.999 percent pure. Resistivity ratios of 1400, 
for annealed specimens at room temperature and 4.2° K, confirmed this 
purity. 

The physical arrangement of the specimen and the dummy specimen 
are shown in Fig. 1. They are supported in a horizontal position by 
fine 0.0005 -in. -diameter gold wire hangers, two of which serve as poten- 
tial leads. The specimen length defined by the potential leads was 
nominally 2 cm for most specimens. The dummy specimen length was matched 
potent iometrically to the specimen to reduce errors due to ambient tem- 
perature changes in the whole assembly. The success of this matching 
procedure was indicated by the agreement of the quenching results from 
these specimens with other specimens 7 cm long similarly matched. 

During the prequench anneal it was important to keep a reasonably 
uniform temperature over the test section of the specimen. In a series 
of preliminary experiments with the system just described, the tempera- 
ture distribution along the test length was found to remain constant 
within 5° C at the quench temperature. This was determined by both op- 
tical pyrometry at atmospheric pressure and by thermocouple measurements 
at high pressures. Since no serious temperature fluctuations occurred 
along the specimens at high pressures, subsequent quench temperatures 
were measured by resistivity changes of the specimen. The proportion- 
ality between resistance and resistivity (i.e., length to cross section 
ratio) was determined for each specimen by comparison of the specimen 


resistance measured at a known room temperature with the resistivity 
data of Meechan and Eggleston for that temperature. The specimen tem- 
perature was continuously monitored during anneals using potentiometric 
strip chart recorders. During the normal prequench anneal, small cyclic 
fluctuations in temperature of approximately ±3° C were observed at 
higher pressures. A reasonable estimate of the total probable error in 
the measurement of the quench temperature is ±5° C. 

Quench rates were measured as a function of pressure by an oscil- 
lographic technique. The specimen temperature was found to be nearly 
proportional to the potential drop across the specimen during most of 
the quench. The time rate of decay of this potential during the quench 
was observed on a calibrated persistent screen oscilloscope. From the 
quenching curves thus obtained, the initial quench rates (average slope 
over the first 100° C temperature drop) and the half-lives (time re- 
quired to fall to one-half of the quench temperature) were recorded. A 
curve of half-life as a function of pressure (shown in Fig. 3) illus- 
trates the marked increase in quenching rate obtained at high pressures. 
It is significant, however, that the variation with pressure above about 
1 kilobar is small. The half-life for the pressure runs in this study 
ranged from 2.5x10“^ to 5xl0“ 3 seconds corresponding to initial rates 
of 10 4 to 5x10^ °C/sec. It has been shown^“ that such a variation in 
quench rate for quenches below 600° C does not markedly affect the total 
quenched-in resistance of gold. 

The measurement procedure followed throughout most of the work re- 


ported herein was as follows. 


Newly mounted specimens were given a high temperature anneal (tem- 
peratures from 800° to 900° C for 1 hr) followed by a gradual cooling to 
room temperature. Resistivity measurements at room temperature after 
subsequent anneals indicated no additional resistivity change after the 
first anneal. The sample holder was then cooled in the pressure vessel 
to the measuring temperature and, if desired, pressure was applied. Po- 
tentiometric measurement of Rg - Rp was made with the balanced bridge 
circuit previously described. Several series of measurements were made 
over a period of time to determine whether the equilibrium measuring 
temperature was reached. Generally about five such consecutive resis- 
tance differences, constant to within ±1 |iohm, established this. The 
specimen alone was then quickly heated to the quench temperature, held 
there for about 3 to 5 sec (sufficient to establish the equilibrium va- 
cancy concentration-*-^), and then quenched by turning off the heating 
current. The value of Rg - Rp was again measured after the system re- 
turned to equilibrium as before, once more with the requirement of con- 
stancy of readings. The change in the resistance difference from the 
value obtained before the quench was then interpreted in terms of the 
excess vacancy concentration characteristic of the particular quench 
temperature and pressure used. 

RESULTS 

To determine the validity of the experimental techniques used, the 
additional defect resistance quenched into several specimens at atmos- 
pheric pressure was determined as a function of quenching temperatures. 
From Eq. (3) it follows that from the slope of a plot of log (AR/R) 
against l/T, the formation energy for the defect may be obtained. Fig- 
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ure 4 shows the fractional resistance quenched into three specimens in 
the temperature range from 510° to 890° C. The curves shown are dis- 
placed vertically for the sake of clarity. The curves for spec im ens A 
and B actually involved the use of the same wire hut with a slightly 
different heat treatment prior to quench. The data from these curves, 
as well as those for specimen C, actually superimpose, indicating the 
effectiveness of the annealing procedures. Some curvature is seen in 
the temperature region above 725° C, indicating that annealing is occur- 
ring during quenches from above this temperature. The observed linearity 
below this temperature supports the assumption that simple defect produc- 
tion predominates. The best fit of the data below 725° C gives formation 
energies of 0.90, 0.93, and 0.97 ev. The average of these three determi- 
nations is 0.93 ev. 

The results obtained from quenches made under high pressures are 
given in Table I. The quenched- in resistance was corrected to 600° C by 
use of an Arrhenius type equation using the value of 0.93 ev for the 
formation energy and the assumption that the change of the derivative of 
the energy with pressure is temperature independent over the small tem- 
perature interval in which the data were obtained. 

A semilogarithmic plot of the corrected resistivities as functions 
of pressure is shown in Fig. 5. The slope of the curve of Fig. 5 in 
Eq. (4) yields a formation volume of AV^ =5.4510.7 cm^/mole. The 
molar volume of gold at 600° C can be estimated from the molar volume 
at room temperature, the thermal expansion, and the compressibility. 

The value obtained is 10.45 cm^/mole. The specific volume change per 

5.45 cm’Vmole 

10. 45 cm^/mole 


defect is then 


= 0.5210.07 at* vol. 
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DISCUSSION 
Formation Energy 

The energy for defect formation obtained from Fig. 4 is in good 
agreement with the values, quite generally accepted as the single va- 
cancy formation energy, obtained from atmospheric quenching experiments 
and from calorimetric, dilametric, and other methods. The best esti- 
mates from these measurements range from 0.94 to 0.98 ev. Above about 
725° C, deviations from linearity appear in the curves of Fig. 4 possibly 
indicating annealing during the quenches from above this temperature. 

This deviation is consistent with the quenching studies of Kauffman and 
others-*- for the quench rates obtainable at atmospheric pressure in this 
experiment. The agreement of the formation energy calculated from the 
lower temperature data with the existing best estimates engenders confi- 
dence in these techniques. This confidence is extended to the results 
from the high pressure quenches, performed at variable but higher quench 
rates. 

Formation Volume 

The value obtained for the formation volume is in good agreement 
with the value of 0.53 at. vol. obtained by Huebener and Homan® by high- 
pressure quenching. It is also in agreement with the value of 0.57 at. 
vol. obtained by DeSor bo, -*--*- who derived AV^ indirectly by combination 
of the vacancy formation energy and his value of the vacancy concentra- 
tion obtained from the energy released on annealing vacancies in gold 
foils with the resistivity change of gold wires per fractional volume 
change measured by Bauerle and Koehler. -*- < - ) 
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Calculated estimates of Tewordt for the vacancy formation volume 
in copper ranging from 0.47 to 0.55 at. vol. 21 are also in agreement 
with the measured value. 

A theoretical estimate of the degree of relaxation of nearest 
neighbor atoms to the vacancy can he made assuming the Eshelby elastic 
continuum approximation. 22 For f.c.c. metals, the relaxation 5 is 
given by 



where Vf is the measured formation volume and a is Poisson's ratio. 

Poisson's ratio of 0.42 for gold 2 ® and' the measured formation volume 
of 0.52 at. vol., give the inward relaxation of nearest neighbors to the 
vacancy as 2.4 percent. This value compares favorably with a relaxation 
of 2.24 percent for Cu calculated using a Morse potential function,, 12 .. 
This remarkable agreement tends to justify the validity of the use of a 
simple pairwise potential function for lattice calculations in which 
drastic atomic rearrangement does not occur. 

The measured formation volume, when added to Emrick's vacancy mi- 
gration volume for gold® yields an activation volume for self -diffusion 
in gold of 0.67 at. vol. This value may be compared with 0.7 at. vol. 
measured for self -diffusion in lead® and 0.9 at. vol. reported for self- 
dif fusion in silver.® To the author's knowledge, the activation volume 
for self -diffusion in gold has not been measured; however, because both 
the formation volume and the migration volume have been measured, a di- 
rect comparison with the diffusion volume would be most interesting. 
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Finally, the Bauerle and Koehler resistivity change of gold wires 
accompanying a fractional volume change Ap/(AV/v) = 3. 2x10" ^ ohm-cm 
and the fractional formation volume measured can he used to calculate 
the resistivity per atomic percent vacancies. The value obtained, 

1.67 ohm-cm/at. percent vacancies, is in fair to good agreement with 
several theoretical estimates for copper made since 1953 24 ranging from 
1.25 to the most recent value 1.67 ohm-cm/at. percent vacancies. 
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TABLE I. - DATA FOR PRESSURE RUNS [AR/R, FRACTIONAL 
RESISTANCE QUENCHED INTO SEVERAL SPECIMENS AT 
TEMPERATURE T Q ; (£R/R) corr , RATIOS 
CORRECTED TO 600° C 


Run 

P, 

bars 

t q , °c 

(AR/R)xl0 4 

(AR/R) corr xl0 4 

1 

0 

509 

1.60 

8.53 

2 

783 

592 

7.35 

8.16 

3 

2300 

612 

8.88 

7.50 

4 

2550 

590 

6.17 

7.18 

5 

3050 

598 

6.68 

6.88 

6 

3720 

570 

4.24 

6.60 

7 

3800 

584 

5.17 

6.48 

8 

4600 

579 

4.42 

6.15 

9 

5360 

591 

5.44 

6.03 

10 

6100 

581 

4.22 

5.60 










E-2413 


! 


Low thermal connector— ^ 



Figure 1. - Exploded drawing of specimen enclosure and cold thimble showing 
arrangement of specimen and dummy specimen wires. 



Figure 2. - Schematic diagram of circuit used for resistance measurements. 
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Figure 3. - Influence of pressure on time required for temperature of speci- 
men wires to fall to one-half quenching temperature, 600° C (ti/2>. 
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